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ABSTRACT 

Time evolution of unstable particles that occur in the expanding universe is 
investigated. The off-shell effect not included in the Boltzmann-like equation is 
important for the decay process when the temperature becomes much below the 
mass of unstable particle. When the off-shell effect is taken into account, the thermal 
abundance of unstable particles at low temperatures has a power law behavior of 
temperature T, jji-^)"^^ unlike the Boltzmann suppressed e~*'^/^, with the power 
a related to the spectral rise near the threshold of the decay and with T the decay 
rate. Moreover, the relaxation time towards the thermal value is not governed by 
the exponential law; instead, it is the power law of time. The evolution equation for 
the occupation number and the number density of the unstable particle is derived, 
when both of these effects, along with the cosmic expansion, are included. We also 
critically examine how the scattering off thermal particles may affect the off-shell 
effect to the unstable particle. As an application showing the importance of the 
off-shell effect we compute the time evolution of the baryon asymmetry generated 
by the heavy X boson decay. It is shown that the out-of equilibrium kinematics 
previously discussed is considerably changed. 
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1 Introduction 



There are many short-lived particles that have existed in abundance in the early 
universe whose temporary presence did not leave behind any measurable effect. Im- 
portant exceptions to this exist, such as the neutron which certainly is the key for 
explanation of element abundance of the present universe. 

Theoretical estimate of the abundance of these unstable particles after the cos- 
mic temperature drops below the mass of the unstable particle is very important for 
subsequent time evolution. Most works in the past are based on the Boltzmann 
equation that takes into account relevant reactions in the expanding universe. The 
use of the Boltzmann equation has however been questioned recently ; a more pre- 
cise quantum mechanical description of the decay process in thermal medium should 
contain important off-shell contribution not properly treated in the Boltzmann ap- 
proach. These off-shell effects are eminent in the low temperature region. Low 
temperature effects are clearly important in this problem, since unstable particles 
are typically very non-relativistic when they disappear in the early universe. 

In the present work we shall develop a general formalism of computing time 
evolution of the net number density of unstable particles and clarify the off-shell 
effect. The off-shell effect appears in two ways; first, in a slower relaxation towards 
the equilibrium abundance and second, in a larger equilibrium value not suppressed 
by the Boltzmann factor such as ^-^^^/t -\vhere AM is the mass difference of 
the parent and the daughter particles. It is shown below that the off-shell effect 
becomes dominant below some temperature Teq. The abundance of unstable particles 
then follows the power law; ^ ~ where a is a parameter related to 

the threshold behavior of the spectral function for the decay and F is the decay 
rate. Thus, unstable particles do not disappear suddenly. Instead, their abundance 
gradually decreases with a power of decreasing temperature as the universe expands. 
Physical processes that follow after the decay are then prolonged. The off-shell effect 
turns out to be more prominent for a larger decay rate. 

We next consider as an illustrative application of this general result the hypo- 
thetical X boson decay that may have created the matter-antimatter asymmetry 
when they decay 0, We find that the time evolution of the baryon asymmetry 
is substantially changed and the severe lower bound of the X boson mass is consid- 
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erably relaxed by the off-shell effect. For the first time we find that some mode of 
the X boson decay for baryogenesis is excluded due to the off-shell effect. This is the 
S-wave decay mode into a boson-pair. 

This paper is organized as follows. In Section 2 the theoretical model of unstable 
particle decay is explained. This is a field theoretical extention of the harmonic 
model for the quantum dissipation in thermal medium discussed in 0. We first 
present and formally solve the quantum mechanical model of the decay of excited 
levels in thermal medium. A great virture of this model is that its integrability leads 
to explicit formulas for many quantities of interest. One can clearly see how the 
off-shell effect arises in these formulas. Extention to the unstable particle decay in 
field theory models can be made, but it is in general complicated and not readily 
solvable. But fortunately, in a thermal medium far away from the degeneracy limit 
which is relevant in the early universe the decay process is approximately described 
by this class of solvable quantum mechanical models extended to infinitely many 
decay channels. In Section 3 the occupation number and the number density of a 
species of unstable particles is calculated and its time evolution equation is derived 
in the expanding universe. The stationary abundance when the cosmic expansion 
is switched off is worked out, and its behavior at both high and low temperatures 
is studied in detail. In Section 4 we pay a special attention to the off-shell effect 
and its role in cosmology. We also discuss a possible effect of the incoherence due 
to the scattering off thermal particles and its role to the decay process in thermal 
medium. In Section 5 we apply previous results to the problem of baryogenesis. 
The time evolution equation for the baryon asymmetry is derived, including the 
off-shell effect. This equation is analyzed both analytically and numerically, and a 
comparison is made when only the on-shell contribution is retained. 



2 Model of unstable particle decay 



We first present an exactly solvable model of the decay of excited levels in ther- 
mal medium, an extention being made to the case of many decay channels. This 
is a slight extention of our previous harmonic model 0. We then explain how the 
two-particle decay of unstable particles in any quantum field theory can be approx- 
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imately described by this class of quantum mechanical models of infinitely many 
decay channels. The approximation is valid for the thermal medium of the low 
occupation number, the circumstance far away from the degeneracy limit. 

We assume that an excited state |1) of energy Ei is given by applying a creation 
operator to the vacuum |0); |1) = c^O) . There are continuously many states de- 
generate with this, 6"'^(ti;)|0) as decay states. Here u is the energy of the continuous 
many levels. Without specifying the nature of the decay process, one may take the 
Hamiltonian that governs the decay as 



H = Eic^c+ r dujujh\uj)h{uj)+ H duj^a{uj) [h\uj)c + c^h{uj)) . (2.1) 



Here ujc is the threshold of the continuous states taken to be Uc < Ei , and o"(a;) 
characterizes the decay interaction. This Hamiltonian is a general one with regard 
to the decay process, complications being hidden in identification of the composite 
operator b'^lu) and the spectral form of interaction a{uj). 

As emphasized elsewhere 0, the dynamical system thus specified is exactly solv- 
able; one may explicitly construct the diagonal operator B''{uj) and the eigenstate 
|ci;)5 = i?"''(ti;)|0) that diagonalizes the decay Hamiltonian; 



a(u)a(u') 



B\uj) = b\uj)+F{uj + tO+) -Ja{uj)c^+ / du b\uj') . (2.2) 

Here F{z) is analytic except on the branch cut along the real axis z > uj^ and is 
given by 

F{z) = i ^ . (2.3) 

J- JUJc LU—Z 

One can explicitly check that the canonical commutation relation, 

[B{uj),B^{uj')]=S{iu-uj'), (2.4) 
and the important inversion relation and the Hamiltonian equivalence. 



b\uj) = B\u) + / du' ^ ' ' ' \ \ B\J) , (2.5) 

POO I 

c^ = - duj^a{uj)F*{u + iO'^)B\Lo), (2.6) 

POO 

H= / du;u;B^Lj)B{Lj). (2.7) 
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The basic reason of integrability is saturation of the unitarity relation by "elastic" 
one; 

F{u + iO+) - F{uj - iO+) = 2m a{uj)\F{uj + iO+) ^ = 2ni H{u) . (2.8) 

The quantity H{u) is characterized as the overlap between the prepared state c^\0) 
and the eigenstate |<^)5; 

H{uj) = mc\uj)s\'. (2.9) 

In the weak coupling of a{u;) <^ M, the spectral function H{u;) has a Breit-Wigner 
form as seen from the formula, 

where n(a;) is real and given by the dispersion integral, 

U{z) =V duj^^. (2.11) 

As is well known, there is a simple pole of F{z) in the second Riemann sheet 
near the real axis which describes the time evolution in the form of the exponential 
decay; e"*^^*""'"*/^ . The imaginary part of this pole coincides, in the weak couphng 
limit, with the decay rate given by perturbation theory; F = 27ra{Ei) . 

Since the operator solution is known, one can explicitly write down many quan- 
tities of interest; for instance the non-decay amplitude of a pure unstable state c^|0) 
is given by 

(l|e-*^*|l) = (0|ce-'^*ct|0) = / duo a{uj) \F{uj + tO+)\^ e-'""' . (2.12) 



We are primarily interested in the occupation number and its time evolution in the 
cosmological thermal medium. Let us first recall this quantity in the pure state 
\ip{t))^ e-'"* |1) at time t, 

{iP{t)\c'^c\i){t)) = (l|e'^*ctce-'^*|l) = {l\c\t)c{t)\l) . (2.13) 

Here c{t) — e'^* c e~*^* is the Heisenberg operator and in this model 

c\t) = - du^a{u)F*{u + iO+)e''''B\uj). (2.14) 

JllIc 

More conveniently, it is written in terms of the parent and daughter operators, 

roo I 

ct(i) = g{t)c' + i / dio^a{u)h{u,t)e'''H\u), (2.15) 

fOD POD 

g{t) = / du;H{u;)e^*^ / du;\{l\uj)s\'^ e''^', (2.16) 
h{uj,t) = iF*{uj + iQ-^)-ik{uj,t), k(u; ,t) ^ ie-''^* g(t) , (2.17) 
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where a condition k{uj , oo) = is imposed such that the asymptotic value is 



h{uj,oo) =iF*{uj + iO+). (2.18) 

The occupation number in the pure state is thus given by 

{l\c\t)cm) = \g{t)\\ (2.19) 

There is a very useful way to compute the basic function g{t). One can use the 
analytic property of F{z) to express this function as a sum of two contour integrals 
as shown in Fig.l; the first one encircling the pole in the second sheet (Co) and the 
second continuous integral along a complex path Ci, 

9{t) = ^ ) dzF{z)e'^' ^ go{t) + gr{t) . (2.20) 

Physically, this second continuous contribution gi{t) gives the off-shell effect, while 
the pole contribution go{t) essentially gives the on-shell effect. Both terms decrease as 
t —>■ oo, but the Ci integral has a power dependence in contrast to the exponential 
form of the pole term; gi{t) oc , where a is related to the threshold behavior 

of the spectral function, <j{uj) ^ c{u! — Uc)" , near u = Uc- 

There is an equivalent and more intuitive way to separate the on-shell and the 
off-shell contributions. For this we go back to the u integral along the real axis. 



9it) = r' dco , ^ , J^^l , , , e^-* . (2.21) 



Assuming the weak coupling, 11 (cj) <^ Ei , one separates the region of integration 
into the two parts, one around the pole, ~ i?i + Il{Ei) , and the rest of the region 
which is dominated near the threshold, uj = uJc, for a large time. This approximation 
gives 

1 r^c 

g{t) ^ exp[iEit - TTa{Ei)t] + duj a{uj) e"^\ (2.22) 

where Ec is a physical cutoff scale of order Ei. The first term is the on-shell con- 
tribution, while the second is the off-shell contribution. In computation of physical 
quantities in thermal medium there are other energy scales such as the temperature 
T, which may replace this cutoff by the factor e'^^"^. The important message is that 
the off-shell effect at late times is determined by the u integration near the threshold 
region. 
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Extention to the many channel problem is straightforward. We denote the chan- 
nel by an index i and write the decay interaction as 

r du ^ (bl{u;)c + c%{u;)) . (2.23) 

Only the key formulas are quoted; 

+F((^+zo+) ( - v^ct + r dj y: biiuj') ] , (2.24) 

\ k — — / 

^(^) = , (2.25) 

c^{t) = g{t)c^ + t r dujh{uj,t)e''^' \f^i^bi(^)^ (2-26) 
h\{uj , t) = e'^* &1 e-^^* = e^-* ( V^^^(^ ' + ^'l(^) ) 

+ / dJ ( ih{u , t)e-* - , t)e-'* ) , (2.27) 

duj ^ a,(a;) |F(u; + zO+) e^-* , (2.28) 

/i(cj , t) = i ( F* (a; + iO+) - A;(tu , t) ) , (2.29) 

kiu.i)^— [ dz^^e'^'-^^K (2.30) 
27ri JCo+Ci z — a; 

Note that both the analytic function F[z) and the basic g{t) are determined by the 
total strength of the spectral function, 

a{u) = ■ (2-31) 

i 

We now consider the field theory model of unstable particles, denoting the parent 
particle by c and two-body daughter particles by hi 62. The decay interaction of the 
unstable particle of a momentum q is given by the Hamiltonian; 

(27r)6 ^ ^ 1 2;^g^^^-^ 

(2.32) 



- / (^^)'^^(' - - ( *l(*04(fe)c(,-) + (h.o.) ) , 



with g some coupling constant. We may then identify the decay product operator 
and the spectral function as 



y'^ii'^ ^Q)bl{'^ ,q) 



/ ^-^^ (2^)'^'(?"*- ^1 - ^2) Siuj -u,- u,) ^^g== h\{k,)h\{k,) , (2.33) 
o-(cj , g) = 27r ^ 0-^(0;, g) = 

i 

^ (27r)3 8uJiUJ2UJg 

where = 1 for fermions and Fi = for bosons. The decay amphtude Ai should 
be given separately in specific decay models. This can be thought of an extention 
of the discretely many (and finite) decay channel problem to the continuously many 
(and infinite) channel problem. 

The field theory model of unstable particle decay is not exactly solvable, because 
the commutator among the decay product operator, 

[hh,b\bl] = l±{b\h + blb2) (2.35) 

(± referring to the boson or fermion pair), is not the cannonical one [&,&^] = 
1 as in the quantum mechanical model. Thus, the composite operators bl{uj,q) 
introduced by eg. ( p. 331) do not obey the canonical commutation relation so crucial 
to the integrability. But the important case in which the bilinear term in the right 
hand side of the commutator can be neglected and the replacement is made, 

[bib2,b\bl] ^ 1, (2.36) 

is identical to the solvable model in quantum mechanics. This occurs in the cir- 
cumstance under which the thermal medium is very far away from the degeneracy 
limit. In this case the occupation number fi which is the expectation value of the 
operators bjbi in the thermal medium is very small (/« <^ 1), and one may neglect 
the bilinear term in the commutator above. Even in the dense, hot early universe of 
the standard cosmology the low occupation number is realized. We shall thus fully 
exploit this approximation in application to cosmology. 

From the definition of the spectral function a relation to the decay rate follows; 

M r 

a(a;,,g-) = — — , (2.37) 

on the mass shell with uj = ujq = \JM'^ + q^. The factor ^ represents the time 
dilatation effect. Off the mass shell, 

M 



a{LU,q) = —criJco^-cp), (2.38) 
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where the spectral function in the right hand side (j{uj) is the one in the rest frame. 

A choice of the decay model corresponds to a particular form of the spectral 
function. For instance, in the fermion-pair (of equal mass m) decay of a scalar boson 
if ipip described by a Lagrangian density of C = g ipipip , 

a{ij , k) = {uo^ -k''){l- ^' ■ (2-39) 

A more general, convenient parametrization of the spectral function that becomes 
adequate in the temperature range of T 2m (the threshold for the decay product 
pair) is given, using the decay rate F; 

a(uj,k) = — ^ — . 2.40 

For instance, the gauge X boson decay into a fermion-pair as well as the scalar X 
boson decay given by ( |2.39| ) has this form with a = 2, while the Higgs X boson decay 
into a boson-pair has this form with a = 0. 

In reality, the spectral function in any non-trivial field theory model is com- 
plicated beyond the lowest order of perturbation. But it turns out that what is 
important for our subsequent analysis is the on-shell value of the spectral function 
given by the decay rate F and its behavior near the decay threshold, hence the pa- 
rameter a in addition to the decay rate. The value of a is dictated by the unitarity 
relation for the opening channel, thus is essentially of kinematical origin. Both inter- 
mediate and late time behaviors of the decay, and important temperature dependent 
off-shell effects are described by these two parameters. Hence in a sense the detailed 
specification of a field theory for the decay is unnecessary. 

The occupation number of the unstable particle decaying in vacuum is thus given 

by \giq,tW with 

^(g,t) = ^== / d^a(v/a7^)|F(^ + ^0+,g-)re^'^S (2.41) 

where are masses of decay product particles. 

We should mention a limitation of our approach. We neglected in our model 
Hamiltonian the interaction of unstable particles with the medium, except the decay 
interaction. This guarantees a coherence of the decay and its inverse interaction 
and makes it easy to treat the environment effect on the decay process. There 
may however be an important class of thermal interactions on the unstable particle; 
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the scattering off thermal hght particles. If the thermal interaction of this sort is 
included, it gives rise to an additional term to the spectral function. In terms of the 
coupling strength this is a higher order effect of order o? T compared to the decay 
rate of order aM. However, this might contribute in the off-shell region. We shall go 
back to this effect when we discuss the off-shell effect in thermal medium in Section 
4. 

In the field theory of unstable particle decay one needs to perform renormaliza- 
tion. The method of renormalization is explained elsewhere |^ and here we shall 
write for simplicity quantities without renormalization, since in lowest order of per- 
turbation renormalization is straightforward. 



3 Generalized Boltzmann equation 

In order to discuss the decay process of unstable particles that occur in the cosmic 
thermal medium, one must take into account the presence of medium and incorporate 
the inverse process that creates the unstable particle. A thermal environment is 
described by the density matrix denoted here by pi and one has to consider how the 
decay proceeds in this mixed state. For the time being we assume that the change of 
the thermal environment is minor and the back reaction of the environment change 
against the decay and its inverse process is negligible. Later when we extend our 
analysis to baryogenesis, we incorporate a relevant effect of the environment change. 

In the mixed state the occupation number of an excited level is 

/(t)=tr (ct(t)c(t)p.) . (3.1) 

It is not difficult to show from the operator solution that this quantity obeys the 
first order differential equation; 

at g 



+ 



duj / du' ^a{uj)a{uj') ( gh* {J , t)e~^"* + g*h{uj , t)e*'^'* 



+ iioj - uj')h(uj , t)h*(ij' , t) - 23fJ^ h(uj , t)h*(u' , t) 1 e'(— ')* (bUiuWuj')), . (3.2) 

9 
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We used the notation, 

(A)i = tr(Ap,). (3.3) 

Although one can write f{t) in an integrated form, this differential equation is more 
useful when one incorporates effect of the cosmic expansion. Another advantage of 
this form is that the initial state dependence via (c^c)j is eliminated in favor of the 
occupation number f{t) at any time t. It is however convenient not to eliminate 
the initial state dependence of environment variables, (6^(a;)6(a;'))j , when we later 
incorporate the environment change. 

We take the uncorrelated initial state satisfying 

(etc), = /(O) , {c^b{uj)), = , {b\u;)b{u;% = fM 5{uj - J) . (3.4) 

Although the unstable particle may or may not be in thermal equilibrium at much 

higher temperatures in the earlier epoch, the choice of the uncorrelated initial state 
seems reasonable, because there exists a time lag between the unstable and the decay 
product particles due to different interaction among themselves and with the rest 
of the environment (not written in the Hamiltonian above). We imagine that when 
unstable particles arc about to decay, the decay interaction cannot keep pace with 
fast thermal interaction in the bulk of medium, hence we assume that the unstable 
particle has no thermal contact with the medium when we start calculation of the 
abundance evolution. In the end of the next section we however estimate how the 
scattering off thermal particles may affect the decay law in thermal medium. 
Introducing the rate defined by 

T{t)^-2^^ = -^ln\g{t)\' (3.5) 
g dt 

(which reduces to a constant rate F in the pole dominance approximation), one has 

r duj a{uj) ( 2^{g{t)h*{uj , t)e-"^') + r{t) |/i(cu ,t)\^) fi{uj) . (3.6) 

It is instructive to first give the time evolution in the narrow width approximation 
and then explain an improved approximation incorporating the Ci integral. Using 

go{t) ^ e-^*/^+^^i* , (3.7) 
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one has 



■ ( r( 1 - e~^*/2 cos(u; - Ei)t ) + 2e-^'/\uj - E{) sm(u; - Ei)t ) . (3.^ 



There are oscillatory terms of frequency of ~ l/(c<j — i^i), but they are averaged out 
by the uj integration. From eq. ( |3.6| ) the fundamental equation in thermal medium is 
found to be 

f = -r(/-/.(i?i)), (3.9) 

in the narrow width approximation. Thus, this equation describes the relaxation 
towards the environment value by the constant rate F. 

The evolution equation (|3.6|) correctly describes the relaxation process of unstable 
particle in medium beyond the narrow width approximation. Since g{t) and 
h{uj ,t) approaches an equilibrium value as t ^ oo, the stationary value of the 
occupation number is 

/■oo 

f{t) ^ / dLua{Lu)\F{LU + tO+)\^f,{iu) = f^. (3.10) 

If one takes a thermal distribution for the initial density matrix, 

f°° cr(u;) 1 

^"^^L (a; - El + n{u;)y + {na{Lu)y ' ^^'^^^ 

The narrow width approximation cr(ci;) ^ i^i of this gives the well-known result; 

but this form is only approximate and is not a good one in low temperatures. In low 
temperatures of T <^ Ei , the uj integral for f^o is dominated by the contribution 
near the threshold; taking the form of the threshold rise cr{uj) ~ c (cu — tUc)" gives 

^-(^/:-^- '-^) 

When T ^ or T ^ ujc, this further simplifies to 

cF(a + 1 

where F(x) is the Euler's gamma function, and 



/co« ' e-^^-T"*', forT«i^,, (3.14) 



cC(a+l)r(a+l) ^„,._ forB,»T»... (3,15) 

{El - uj^Y 
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with ({x) the Riemann's zeta function. The Bohzmann suppressed temperature 
dependence of the occupation number e~'^^/^ near T = Ei is thus changed to the 
power behaved T"+^ at low temperatures 0. What caused this big change is that 
the full Breit-Wigner shape is cut off effectively at the temperature T, since T -C 
the center location of the Breit-Wigner function. 

We now turn to the unstable particle decay. The occupation number for a mode 

k is 



f{k,t) = \gik,t)\'UiLUk)+ J duj\hiuj,k,t)\' 
I ^^M^^^iH^ ^ \M\'U{k,)U{k,)5\k-k,-k,)5{u-h-k,). 

J (27r)3 %kik2UJk 

(3.16) 

— * 

We took the massless particle for decay products such that uji = \ki\. The stationary 
limit of the occupation number is then 

d?kiCp'k2 



(27r)3 

E \M\'U{h)ft^{k2)5\k-h-k2)5{uj-k,-k2). (3.17) 



foo{k) = I duj \h{uj ,k ,oo)\^ 

g-" n ^ (2cu,)^' 
8/ci/c2a;fc 

In the low occupation number limit this is further simplified since 

/th(A;i)/th(fc2) = e-^(^^+'=^), 
which is replaced by e~^^ in the above integrand. It is thus found that 

/.(A?) . f -^-^^^)--'^ _ (3.18) 

A (cu-cu,,)2 + 7rV2(cu,A;) 

Separating the pole and the threshold regions for this lo integral, one has 

foo{k) ^ e-"'^^ + ^-^j-, / dx {x" - ir'^ e-"''^ . (3.19) 



27rM"+2 ji 

The last integral in this equation is given by the modified Bessel function, and in 
the large k/T limit it behaves as 

r/^e-'^/^. (3.20) 



This has a narrower spread of the momentum than of order \fMT for the Maxwell- 
Boltzmann distribution, and its average is of order T (actually (| + 3) x T using the 
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distribution, ( p.20| )). We shall later use this fact to simplify the momentum depen- 
dence of the effective decay rate at late times. We plot the momentum distribution 
given by eq.( |3.19| ) in Fig. 2 to compare with the on-shell distribution relevant at low 
temperatures; e~ V'^^^+a^^/^. 

The total number density n{t) is then obtained by summing the occupation num- 
ber over modes of various particle momentum k. In an isotropic medium it is 

1 r'^ 

n(t) = ^ i) ^^^'•/'(^'^)- (3-21) 



In general, various quantities in eq.( |3.6| ) depends on k. This dependence is traced 
to the boost of the parent unstable particle. We shall often suppress the k depen- 
dence, unless otherwise there is a confusion. The stationary number density rioo is 
determined using f^o for f{k in the two-particle decay of equal mass m 

n^ = ^ r dke r duj , ^^"" '^^ .(3.22) 

Throughout this work we ignored a small part of n(ti;); the residual n(ci;) that 
remains after renormalization. 

In the low temperature range of M ^ T ^ 2m the time dilatation effect is 
negligible and 



oo 



^ , , duj—. / dkk'^ — -— ^ — . (3.23) 

47r3 Jo e^^ -1 Jo M"+2 ^ ^ 

It is then analytically calculated as 

^ ~ A(a) ^ i^r^' Ma) - C(« + 4)r(a + 4)r(f + i) 

In Fig. 3 we plotted the stationary number density, eq. ( |3.22| ) combined with the 
spectral function ( |2.40| ), as a function of T/M for two values of the decay rate F. 
Our numerical computation supports that to a good accuracy the stationary number 
density rioo is given by a sum of the pole and the threshold contribution; 

1 A"" A;2 ^, V ,T 



This formula is accurate for any temperature T less than M if the decay rate F/M is 
small enough. But at higher temperatures the first on-shell term in eq. ( |3.25| ) alone is 
accurate and the off-shell power term oc T"^^ in this equation should be discarded. 



14 



One may define the equal temperature Teg as the one at which the two contribu- 
tions in this equation, the pole and the threshold contributions, become equal. How 
this equal temperature depends on the rate F/M is important, and it is shown in 
Fig. 4 for two values of a. With the on-shell contribution given by 

noo^(^Pe-^/^ (3.26) 
the equation that determines Teq is 

a;"+ie-^ = (2vr)='/M(«)I:, x^^. (3.27) 
for a small F/M. A rough analytic estimate in the F ^ limit would be 

but this expression is accurate only for a very small for instance, j-^ < 10~^ for 
a = 0. 

We now turn to the time evolution of the occupation number and the number 
density. The approach towards the stationary value is governed by the large time 
limit of g{t). In the pole dominance approximation the relaxation is exponential in 
time and fast. On the other hand, the true late time behavior is the power law, and 
this makes complete relaxation slower. p[ 

The effect of the cosmological expansion is readily incorporated for the number 
density n by adding the 3H n term where the Hubble rate 

a 1 , ^ 

H = - = —. 3.29 

a 2t ^ ^ 

We first give the result in the narrow width approximation; 

^ + 3-n = -r(n-n*^(T)), (3.30) 
at a 

n^^{T) = — dx , : . (3.31) 

We introduced an averaged quantity over the momentum in order to simplify the 
mode integral; 

1 roo 

r = 7r^7TV / dkeT,f{k,t). (3.32) 
Ztt^ n[tj Jo 

In low enough temperatures this averaging is indeed simple; F ^ F , because the 
time dilatation effect given by M/a/M^ + ^ 1 — ^ is negligible for the slow 
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motion of parent particles. The equation for the number density should be combined 
with the well known temperature-time relation in the radiation dominated universe, 



drripi a T 45 

with N the number of massless species contributing to the energy density. 

In the low temperature region of T < M the equation is further simplified by 
using the Maxwell- Boltzmann distribution function of the zero chemical potential, 

g-(M+|i)/r_ (334^ 
A convenient evolution equation is obtained using the dimensionless quantities. 



Y=— u= — 



r^rt. (3.35) 



T \\dmpi' 

The time scale of variation is given by the hfetime F"^, hence it is useful to use the 
dimensionless time r. The evolution equation is then 

^^-(Y-S), S{u) ^ i^f' e-- . (3.36) 



The source term varies as 5" {\J^ i "^i^^ 



.= %^. -,/|-,/^. (3.37) 



Thus, 5" as a function of r has a maximum around t — rj. The meaning of the 
quantity rj is the decay rate F divided by the Hubble rate H at the temperature 

r = M. 

One can readily integrate this equation. The solution to the differential equation 
for the yield Y is 

Y{r) = £ dr' S {J^ e- (^"^'^ + Y{n) e" (^"^^^ . (3.38) 

The first term gives the yield created from the thermal medium. For r ^ 1, namely 
at times much larger than the lifetime, t ^ 1/F, the dominant region of the r' 
integral is r — r' < 1 and the integral gives 

y{r) ^ S U^^ . (3.39) 
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This is in general valid unless r] < rj^r , where the critical value ?7cr = The yield 
Y thus roughly follows the thermal value S{T). 

There is some lesson one can learn on the more general case of the mode indepen- 
dent r(t) = — In p, from this calculation in the pole dominance approximation. 
It is obvious that even for the more general case the general solution to the mode- 
summed form of eq.(|3.6|) extended to the expanding universe is given by 



mH.Wn/'^*'^ + m)^), (3.40) 



dk k'^ / du 



where the source term S{t) is given by 

1 r- 
27r2 T3r(t) io 

• ( 2^igit)h* ico,k, t)e-^-*) + r(t) \hico ,k,t)\') . (3.41) 

The stationary value of this source term S{t) at t — > oo is 



Soo = ^, (3.42) 

where rioo is given by eq. (|3.22| ). Actually, this quantity is not stationary, since the 
temperature T gradually changes with the cosmological expansion. 

When the rate Tk does depend on the momentum k, the solution for Y cannot 
be given in a simple closed form such as eq. ( |3.4CI| ). 



4 OfF-shell effect 

The off-shell effect appears in two ways; first, in the slower time dependence of 
relaxation {\g\ oc t~"~^), and second, in a larger source term (5* oc T""*"^ at small 
temperature T). Both are related to the threshold behavior of the spectral function, 
cr oc (u; — cJc)" . 

We found elsewhere that adding a power term gi(t) to the pole term gQ{t) is 
an excellent approximation for g{t) in the entire time range, unless one looks into 
the very short-time region of relaxation. The power law period is represented by 

9i(t) = ' gl^at/^ e-(--*+-/^) , k,{uj,t) = -zf dt' g,{t')e-^-'' , (4.1) 
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where cr{uj) ~ c{uj — uJc)°' near the threshold and Q = M — uJc ■ With this power 
behavior, the rate r(t) = - | In \g{t)\'^ is ^ 2(a + l)/t . The final yield Y = n/T^ 
is then of order, Soc oc {'m)"~^^ ■ Thus, the yield does not decrease as rapidly as 
might have been expected from the exponential decay law, but it decreases with a 
power, T"+i oc t-("+i)/2 . 

A complication arises when one incorporates dependence on the particle momen- 
tum. This effect appears in two ways; first, in the time dilatation of the lifetime 
uJk/M and second, in the function g{k,t). The time dilatation effect is negligible if 
one only considers the temperature range of T < M (the mass of unstable particle). 
We shall thus discuss the momentum dependence of the off-shell contribution to g; 
gi{k,t). In the Ci contour integral of Fig.l one has an approximate expression of 
the form, with Uc = k. 



S^^^^^'^-'jp i dy<k + iy,k)e~y\ (4.2) 
We then use the spectral form (|2.40|) or its low temerature approximation, 

oiu.k)^—^ ^ — , 4.3 

to get 

,,ik,t)^,.(2,r''n- + i)-j^,:,^. (4.4) 

An approximate evolution equation is summarized using the dimensionless time 
variable, r = Tt, and the time-temperature relation T/M = ^rj/r , 

dY dkP \( it\ f°° J a{uj,k) 1 



dr Jo 2.^T^^^^^'^y'^'^- L (.. - c,)^ '+ (W J ' ^'-'^ 

I dkk"^ ^ , ^ , , 

''"f/o l^rAW, (4.6) 

2 s _4 +,,(,.,)).- . (4.7) 

Thus, the momentum dependence gi{k ,t) is convoluted with other momentum 
dependent functions in the k integral. In low temperatures the most important part 
of this momentum dependence is in the spectral function a{u! ,k) which vanishes at 
the threshold, uj(. = k for the massless daughter particles. In some sample numer- 
ical computations that include the phase factor in gi{k,t), we observe oscillatory 
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behaviors around the transition time from the pole to the power period that occurs 
at ~ Teq given by (|3.28|) . In the rest of the time region the yield Y smoothly varies, 
and it is very well described by using a simplified effective rate, 

\9o{t)r + \9i{k,W 

The oscillatory behavior in the transition region gets smoothed and approaches 7 
given above when one time-averages the rate and increases the resolution time At 
towards the lifetime l/F. Moreover, the late time behavior is insensitive to whether 
one uses the exact rate 'y{k ,t) or the effective rate 'y{k ,t) above. It is thus a 
reasonably good approximation to use the average rate ([4.8|) . 

In a still further simplification one may neglect the momentum dependence in 
\gi{k ,t)\'^ oc A;" and replace the momentum by its average in low temperatures; 
/c° 0[1] X . The relation k°' oc T" is consistent with the late-time momentum 
distribution of the stationary occupation number foo{k), as discussed in eq.( p.20| ). 
We have checked that the following crude rate equation, 

dY 

-j^ = -l{Y-Yo-So), (4.9) 



e-^^ + (a + 2)i?(a)(^)"+Mir(rt)— ^ 
^ e-r* + 5(a) (^)"+Mi)" (H)— ^ ' ^ • ^ 



1 k'^ 



^ C(a + 4)r(a + 4)r(f + l) £ .^..+1 ^ 
° 167r2 0Fr(f + f) M^m' ' ^ ' 



is a reasonably good approximation. We may use -8(0;) = ""l"!^! , since k°' = 
B{a) for the late-time distribution function already discussed. The source term 
is separated into the on-shell contribution Iq plus the off-shell contribution Sq. 

There is a second temperature T* or time at which the off-shell effect becomes 
conspicuous. This is the time when the two terms of g{t) becomes equal; |(7o(^)P = 



\gi{k , t)P . Using the formula above, one gets in the F — > limit 



M ■ 



M y (a + 4) In f 



(4.13) 



Unless rj is small (roughly rj < l/ln(^)), T* > Tgq usually. The transitional period 
then lasts for a while at > T > Tcq . Thus, T* is the temperature at which the 
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slower decrease of the remnant becomes apparant, while Teq is the one at which the 
larger source term becomes visible. 

The late time limit of solution to this equation can be analytically worked out. 
Since the on-shell source term Yq is exponentially small at late times when T < Tgq, 

^«_^(y_S„). (4.14) 

CLT T 

This equation is readily solvable, and in the r — oo limit 

^ a + 3 T3 - a + 3 ^^"^ M " ^^"^^^ 

This asymptotic form becomes relevant, starting at T = Teq. 

We show some numerical results in Fig. 5 and compare with the simple analysis. 
Besides a, important parameters in the time evolution are the rate T/M and rj that 



appears in the time-temperature relation (|3.37|) . It is clearly seen that the yield 
Y = n/T^ is accurately given by the on-shell formula at high temperatures and 
by the power law formula ( |4.15| ) at low temperatures, with a transition region at 
-p: ^ (several x 10)"-^ . 

One may summarize this result, by saying that the time evolution gives the yield 
Y(t), starting from the stationary value S'oo(T) at high temperatures and ending at 
^^^5*00 (T) at low temperatures. The transient temperature is characterized by the 
temperature of order Teq — T* and the yield below this temperature follows 

Y^Ye^i^r^'. (4.16) 

eq 

Importance of the off-shell effect is measured by how close the temperature 
Max(T^, , Teq) for the onset of the off-shell effect is to the temperarure scale for the 
decay. This temperature is estimated as follows. One first defines T^ by F = H{Td), 
which gives T^ = y/fjM. We then discuss two cases separately. First, if r/ > 1, 
Td > M and the decay and its inverse decay frequently occur at temperatures be- 
tween Td and the threshold M, maintaining the thermal abundance Y = Yq. At 
T < M the inverse decay is less frequent than the decay and there is a Boltzmann 
suppression e~^^'^ until the temperature decreases down to T^. (In the case oir] > 1, 
T=K > Teq always.) Thus, the importance of the off-shell effect is measured by how 
large is the quantity, 

(4.17) 



M Y In p 
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if this number is less than unity. It however often happens that 77 ^ 1 is huge and 
^ M > Teq , as it occurs for instance in the top and the weak boson decay. In 

this case a true measure of the off-shell importance is given by how large the value 

of ^ « l/(lnM/r) is. 

On the other hand, if 77 <^ 1, then <^ M, and the decay does not occur until 

T < Tfii much below the Q value. In this case the off-shell importance is determined 

by how large 

is. 

In these two cases of both large and small 7y, a large decay constant of for example, 
;p = 1 — 10~^ , is expected to give a large off-shell effect irrespective of the -q 
value. This is physically reasonable since for a large coupling the narrow width 
approximation is expected to break down. 

We now turn to the scattering effect mentioned in Section 2. We first give a 
crude estimate and later elaborate more quantitatively. The loss of coherence due 
to the scattering off thermal particles occurs in an inverse time scale of order 

(Js ~ nth(r) 1^ , (4.19) 

where is the averaged cross section of the scattering interaction with the thermal 
medium. This gives a new term to the spectral a of order 

4(2^+1)03) ajT3 

3^2 ^2 ' (4.20) 

where ag is the strongest dimensionless coupling with the medium and 2s -|- 1 is the 
spin degrees of freedom. We took the Thomson type cross section oc a^/M^. This 
is a small addition to the on-shell rate due to the extra coupling factor and the 
temperature suppression. More important is its possible contribution to the late- 
time and the low-temperature region. Its contribution to the occupation number in 
the stationary limit, 

1 

^foo^j^ I duase'^-, (4.21) 



gives 



^ 4(2. + 1)C(3) T 
<^/oo «.(^) . (4.22) 
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We have so far ignored the momentum dependence of the spectral function 
(T<j(ci;,g). Indeed, the above estimate of the spectral function is valid only near 
the mass shell. The correct formula including the off-shell contribution is 

a.(^,g-) = ^ / ^ e-P' d{uj + k - ^{k + q)^ + M^) . (4.23) 

This gives to the stationary number density an extra off-shell term for the scattering 
effect, 

^^00 = ^ J ^ / diue-^'^asiuj^q) ^ ^ t;S,MTnth(T) . (4.24) 

This formula is valid only when deviation from the thermal distribution is small. 
Compared to the on-shell value nth, this is smaller by a factor -p. Thus, 
when the off-shell contribution of order, (p^)"^^ , dominates over the on-shell 
value rith, the scattering effect to the source term S may be ignored. 
The scattering effect also gives a new late-time contribution; 

I f d^q f , , ^ 2(2s + l) 2 T 

(4.25) 

to be compared to the previous off-shell contribution. 

Note that nth(T) < (-^)""'"^T^ in the temperature region of our interest. Thus, 
the scattering effect appears negligible, because the temperature inequality usually 
obeyed. 



T r 

gives a smaller late-time contribution than eq.( [4.26D . 



,.<0[(^V7)^/^], (4.27) 
M rjaj. M 



5 Species evolution equation and baryogenesis 

It is of great interest if one can predict consequences of the modified abundance 
of unstable particles during their decay that become permanently imprinted in the 
rest of the cosmic evolution. In this respect it is important to clarify a possible 
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change of the thermal environment due to the decay. The mere increase of the total 
number of decay products is hardly recognizable in thermal medium. One should 
examine a more detailed distribution of particle species in the process of unstable 
particle decay. As an important example of this class, we shall discuss baryogenesis 
in a simplified model. 

To this end we need to introduce several decay modes distinguished by a flavor 
index j of the continuous state 6^(ti;)|0) . In the problem of baryogenesis we think 
of many channels of different baryon numbers such as the diquark (qq) and the 
quark- lepton pair (ql), considering the decay of anti-X boson along with the X- 
boson. It is presumably better to consider channels of different B — L such as 
N ^ I H , I H (where iV is a heavy right-handed Majorana neutrino, and H is the 
Higgs doublet), in view of that the baryon and the lepton numbers are redistributed 
at lower temperatures by the baryon non- conserving electroweak process keeping 



The basic assumption taken in estimating the environment change is that there 
exists baryon conserving strong interaction among environment particles such that 
the kinetic equilibrium among them is readily established, leading to the environment 
particle distribution described by a thermal distribution function of a finite chemical 
potential /i associated with the baryon number. 



Here aj is the baryon number of the species j. In subsequent application to baryo- 
genesis the limit of a small chemical potential is relevant. Thus, 



5 - L unchanged. §, 0, §, i 




(5.1) 




(5.2) 



It is useful to note the conservation of particle number. 



d 



( 



c\t)c{t) + duY. h]{uj , t)hj{uj ,t) = , 



(5.3) 



which holds as an operator identity. It is easy to confirm that if 




(5.4) 
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as t ^ oo, but a net baryon number may be generated, since 

poo 

/ duj J2 aj {b]{u} ,t)bj{uj ,t)) . (5.5) 



di 



One can write down a set of time evolution equation of each species along with the 
one of the parent unstable particle. These equations are not particularly illuminating. 
We shall directly work out the situation relevant to baryogenesis. In this problem 
the environment is gradually changed by the presence of a small chemical potential 
associated with the baryon number. We imagine that the bulk of thermalization 
processes not included in the decay interaction conserves the baryon number, and 
the time variation of the net baryon number is driven by the asymmetry generated by 
the pair decay of X and X. One thus expands in power series of the small chemical 
potential and identify the baryon number density of the thermal environment; 

^ --,T.^1 r dkk'-^-^^. (5.6) 
TT^ j ^ Jo dujk eP^i' - 1 ^ ' 

Furthermore, we only consider for simplicity the case of many decay modes whose 
rates differ only in the overall partial rate, (yj{ijj) = •yjcr^uj), where the species 
independent (7(00) is the total spectral function common to X and X due to the 
CPT theorem. Thus J2j Ij = J2j 7j = 1 • This approximation is excellent if the 
decay products are much lighter than the unstable particle, as it happens in the X 
decay into ordinary fermions. The fundamental baryon asymmetry e when a pair of 
X and X bosons decay arises via a combined effect of baryon non-conservation and 
CP violation [|IU| such that for some partial rates •jj 7^ 7^ . It is given by 



e = Y.aj (7, - 7j) . (5.7) 



We consider both the fermion-pair decay []TT| and the boson-pair decay |T2[ of X 
boson. It is thus convenient to parametrize the spectral function as in (|2.40|) , 

a{uj,k) = - — , 5.8 

where a = 2 for the fermion-pair decay and a = for the boson-pair decay. 
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To lowest order of the asymmetry e and the chemical potential /i, the basic 
evolution equation is 

-e r du,a(io) {V^(t)\h(u,,t)\' + 2-&(g(t)h'(iO,t)e->"'))f\u,) 



luic doj 



+ "2 dcu (rk{t)\h{u; ,t)\'' + 2U{git)h*{u; ,t)e-'^')) a{u;) f\u;)^ , (5.10) 

^ d „a, , , f°° dk k"^ , ^ , . , ^ ^ . 

a Jo ZTT"^ 



r duj{a{uj) (Tk{t)\h{uj ,t)\' + 2^{g{t)h*{uj ,t)e-'''^')) . (5.II) 



dt 

dcu 
Here 

= i^:;^ (5-12) 

is the thermal occupation number for the zero chemical potential and 

^ = I] "i(7i + 7i) • (5.13) 
j 

We first consider the pole dominance approximation. In this approximation a 
simple relation, 

mp = ^(^-^\h\^ + gh*e-'^'^ , (5.14) 
holds. Using the relation of the baryon number density and the chemical potential, 

nB = \Y.^^jl^T\ (5.15) 



3 . 



the relevant evolution equation is 



| + 3^)n5 = er (n+-n*'^)+r5n_--r^nB, (5.16) 
| + 3^)n+ = -r {n^-nf) , (5.17) 
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< = i^fl'e-^'l\ (5.19) 
dkk'^ , „ „ , dkk"^ 



We simplified the equation, using the condition of very low temperature, T <^ M . 
In the numerical computation below, we allow the temperature region, T ^ M . 
Since the bulk of the cosmic medium is in thermal equilibrium, the usual 



temperature-time relation ( |3.33| ) holds. One can then integrate the equation for the 
average number density (|5.17|) , which gives a source term Ys to the rest of equations; 

= (27r)-3/2 e- ^^"rfa;x=^/^e^-v^-r3/4e-v^)+ 7^3/4 e-(---')(^),, (5.22) 



dYe K 



dr 2 



-YoYb + 6Y_ + e Ys , (5.23) 



'-^--y-^lyoYs, (5.24) 



th 



Solutions to this approximation will be compared to a more precise numerical result. 

A comparison with previous works [|T^,|]T^ reveals some differences. In our treat- 
ment two-body processes 6^(ci;)|0) ^ 6^(a;)|0) are not included. (In the terminology 
of [|13] this corresponds to the off-shell two-body contribution and their on-shell two- 
body term is automatically included in our approach, too.) This is a higher order 
effect, hence was neglected in the present work. The major difference absent in 
the past work and included here is however the off-shell effect, as will be discussed 
shortly. 

An improved, but a still simple approximation that incorporates the power law 
period of decay is to add the two contributions in the exponential and the power 
law periods incoherently, ignoring the interference between the pole and the power 
terms. This introduces a time, or temperature dependence of the rate r{t) that 
includes the off-shell effect. 
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Another simplification is necessary for a practical computation of momentum 
integrals containing the momentum dependent effective rate and the time dilatation 
factor. It is not difficult to partially include the momentum dependence in quantities 
such as 

dkk"^ , , a(uj,k) 1 

JO 2tt^ Jk [uj — uJk) + [n a{uj , k))^ ei^"^ — 1 

Indeed we did perform these momentum integrations in our numerical analysis. But 
unless the integro-differential equation for the occupation number, instead of the 
ordinary differential equation for the number density, is directly solved, it is difficult 
to deal with the momentum dependence in the convolution integral containing fx 
and fx in eqs. (|5.11| ). This is because these distributions are unknown before we 
know the complete answer to the problem, which is a formidable task. Thus, we 
made a replacement of the momentum factor by a temperature factor in those terms 
containing fx and fx] 

using the momentum distribution function relevant at late times. 

It turns out that the rest of the detailed momentum dependence integrated nu- 
merically is not crucial; its variation only changes the transient time dependence 
around the temperature T^, and quantities in the rest of the time region are insen- 
sitive to this momentum dependence. We therefore write down a simplified rate 
equation replacing the momentum dependence by its average above. In terms of the 
net number Y± = n±/T^ and the rescaled time r = Tt, 

!^ = -^{Y+-Yo-So) , 

^ = 1 {eY+ + 6Y^-eiYo + Sq) - f (Fi + ^i) Yb)-c{Y, + S2) Yb , 

(5.28) 

,_ e-^^ + (« + 2)i?(c.)(^mir(n)— ^ r(|c. + 3) 

^ e-r* + 5(«)(^)"^-Mi)Mrt)-- ' ""^"^ r(i« + 3)' ^"-'^^ 
_ C(c. + 4)r(a + 4)r(f + i) r t 

167r2V?r(f + 1) M^W ' 
c^'^^^^^^:^-^K.. (5.31) 

s -gfa)- ^(^ + 3) r(f + i)r(a + 3) r t 
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- - fVTf) M^M^ ' ^^-^^^ 

_ 1 2k^ + 1 

The low temperature approximation was not assumed here for /**^(a;), hence 

1 k"^ 
^'^2^^ Jo eV^-+MVT _ 1 • (5-35) 

It can be readily proved by a rescaling argument that both F_ and the baryon 
asymmetry Yb is in direct proportion to the fundamental CP parameter e. We 
assume that o; < 2 as required for any renormalizable decay interaction. 

Some results of numerical integration of the time evolution equation are presented 
in Fig. 6 and Fig.7. The time evolution for a = and a = 2 is evidently different, 
as seen in Fig. 6. Notably, the final Yb vanishes for a = 0. This difference will be 
understandable analytically, as will be discussed shortly. The final amount of the 
baryon to the photon ratio, approximately given by Yb = is shown as a function 
of r\ in Fig.8, the ratio of the decay rate to the Hubble rate at T = M. In all these 
computations of Fig. 6 — 8 the initial X abundance was assumed to be the thermal 
value. In Fig. 9 we show how the time evolution is affected by the initial condition, 
taking grossly different X boson abundance from the thermal value. In Fig. 6 — 9 
we used for the parameters given by the underlying theory, 5=|, ^ = |) c=6, 
chosen to be consistent with the CPT constraint. 

The asymptotic behavior of solutions to this rate equation differs, depending on 
whether a > | or a < |. First, let us assume that the baryon asymmetry Yb 
approaches a constant (including 0) asymptotically as t — > oo. The first and the 
second equation for Y± in the above set then gives the asymptotic solution; 

2q; + 4 noo 

^-^ITTTit. (5.36) 



3(a + 2) SYb 



Y- ^ \ [ J ^^S^. (5.37) 



Both Y± (X ■^i^)"''^^ The asymptotic equation for Yb differs, depending 

on whether a 7^ 2 or a = 2. We shall first discuss the simpler case of o; 7^ 2. The 
equation for the asymmetry in this case is 

flY 

-^^-cS2Yb, (5.38) 
otr 
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where c is given by eq. (|5.31|) . Since 5*2 oc r (2a+i)/2 ^ 

Yb ~ const. exp[— c / dx S2{x)] (5.39) 

J TO 

approaches a constant for a > |. Thus, as t — > oo, approaches a constant from 
above, its rate to the asymptote being of order 0[t~*^^"~^^/^] . When a = 2, the 
asymptotic equation for Yg is more comphcated, but the asymptotic behavior of the 
asymmetry is identical to the case for 2 > a > |. 

On the other hand, for a < | the integral in the exponent of eq.( p.39|) is divergent 
as r — s> oo, hence Yb — > 0. The asymptotic behavior is as follows: Yb = 0[e~^^^^^'^ "] 
with c > for a; < |, and Yb = 0[t^^^] for a; = |. In superrenormalizable models 
of the boson-pair decay a = 0, and this case gives a vanishing asymptotic baryon 
asymmetry, unless bosons in the decay product quickly decay further into ordinary 
quarks and leptons. 

The amount of the generated baryon asymmetry in the large t] limit is of interest, 
because in the usual treatment of the out-of equilibrium condition a very large mass 
is demanded for the X boson due to the on-shell kinematics [Q: it requires that 
the ratio of the decay rate to the Hubble rate at the temperature T = M, precisely 
rj, should not be too large in order to get a sizable X abundance for the baryon 
generation. More precisely, the on-shell Boltzmann approach gives the asymmetry 
that depends on rj like oc r^"^'^ for a large rj [14|. In Fig. 8 we observe that the 
off-shell effect gives a less decreasing behavior as rj increases. It would be of some 
interest if one can work out the infinite rj limit analytically. 

According to the recent analysis of the reheating problem after inflation it is 
quite possible [|l^ that the heavy X boson has been created right after the explosive 
decay |T6| of the inflaton oscillation and prior to the thermalization process. It 
is thus of considerable interest to examine the baryogenesis without assuming the 
thermal abundance of the heavy boson initially. The non-thermal initial condition 
has been taken in Fig. 8 and for a large enough rj the final baryon asymmetry is seen 
insensitive to the initial condition. 



In summary, we derived the cosmological time evolution equation for the abun- 
dance of unstable particles, including the off-shell effect not taken into account in 
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the Boltzmann approach. Apphcation to the baryogenesis problem shows that the 
out-of equihbrium condition based on the on-shell kinematics is changed. 
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Figure caption 



Fig.l 

Contour for the uj integral that separates the on-shell contribution (Co) from the 
off-shell one (Ci). The dashed parts are in the second Riemann sheet continued via 
the cut starting from the threshold uo^. 

Fig.2 

Momentum distribution for a temperature and a decay rate, T = 0.07 M, F = 
0.01 M . A comparison is made between the on-shell distribution function, e~ v^^=^+m^/t 
appropriate at low temperatures, and the off-shell distribution of a = and a = 2. 
Two of these distributions should be multiplied by 0.01 to get the correct values. 

Fig.3 

Stationary number density given by eq.( p.22| ) which is divided by temperature^ 



for two values of the decay rate, = 0.1 , 0.01 . The dotted lines are calculated using 
the approximate formula, eq.( p.25| ), while the broken line is the on-shell contribution 
alone, the first term in this equation. 

Fig.4 

The equal-temperature at which the on-shell and the off-shell contributions be- 
come equal is shown for two values of a = , 2 . The dashed and the dotted lines 
are result of the approximate formula, eq.( |3.28|) . 



Fig.5 

Time evolution of the yield Y = ^ for different values of r^; 0.1 (Fig. 5a) and 1 
(Fig. 5b) and for different decay rates, = 0.1 , 0.01 . For comparison the on-shell 
evolution (the broken line) and the late time form, eq. ( |4.15| ) (the dotted line) are 
shown. 

Fig.6 

Comparison of the time evolving baryon asymmetry. The case of a = shown 
by the solid line and enlarged in the inlet gives the vanishing value for the final 
asymmetry, unlike the a = 2 case shown by the dotted line. 
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Fig.7 

Time evolution of the baryon asymmetry. Two cases of different decay rates, 
^ — 0.1,0.01, are compared to the evolution given by the on-shell contribution 
alone (the broken line). In the inlet detailed behaviors are stressed. 

Fig.8 

Final amount of the baryon asymmetry plotted against r) {— decay rate/ Hubble 
rate at T = M) . For comparison the result based on the on-shell contribution alone 
is shown by the dashed line. Those marked by open boxes and circles are results for 
smaller decay rates 

Fig.9 

Dependence of the asymmetry on the initial X abundance for different rj values; 
0.1 (Fig. 9a) and 10 (Fig.9b). Result for the initially thermal abundance given by the 
solid hne is compared to those of 10 times the thermal value (the dotted line) and 
the zero abundance (the dashed hne). 
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